NASA Contractor Report 3154 



i ■ 

COPYl~RgEE>BW ?^ S 
AFWU TECyriAUC^ba^l^^S 
^ KlRTJ^NO,A0asS-"^’ "" 


Numerical Optimization Techniques 
for Bound Circulation Distribution 



for Minimum Induced Drag of 
Nonplanar Wings: Basic Formulations 


John Kuhlman 


GRANT NSG-1357 
JUNE 1979 





TECH LIBRARY KAFB, NM 



NASA Contractor Report 3154 


TECH LIBRARY KAFB, NM 

I 


□Dtieii 


Numerical Optimization Techniques 
for Bound Circulation Distribution 
for Minimum Induced Drag of 
Nonplanar Wings: Basic Formulations 


John Kuhlman 

Old Dominion University Research Foundation 
Norfolk, Virginia 


Prepared for 

Langley Research Center 

under Grant NSG-1357 


IVIASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Office 


1979 


SUMMARY 


A theoretical method has been developed for determining the optimum span 
load distribution for minimum induced drag for subsonic nonplanar configura- 
tions. The undistorted wing wake is assumed to have piecewise linear variation 
of shed vortex sheet strength, resulting in a quadratic variation of bound cir- 
culation and span load. The optimum loading is obtained either through a 
direct technique, whereby derivatives of the drag expression are calculated 
analytically in terms of the unknown wake vortex sheet strengths, or through 
use of Munk’s criterion. Both techniques agree well with each other and with 
available exact solutions for minimum induced drag. 

INTRODUCTION 

One way in which to improve aircraft performance at subsonic cruise is to 
utilize nonplanar lifting surfaces designed for minimum induced drag. At least 
two theoretical design methods have been developed for this purpose; one by 
Lamar (ref. 1) and another by Feifel (ref. 2). Ishimitsu (ref. 3) uses Feifel’s 
technique for the design of winglets. Both of these theories use vortex lattice 
representations on the wing with trailing filaments extending streamwise into 
the wake. Feifel notes in reference 2 that the discretized vortex lattice 
technique can lead to appreciable errors in local velocity, and hence in the 
span load, for nonplanar configurations in the vicinity of an abrupt change in 
dihedral. Even though the effect on total coefficients is negligible, camber 
solutions in such regions can be in error. One way in which to minimize such 
problems is to assume a more complicated functional form of the bound circula- 
tion than piecewise constant. This avoids the isolated singularities of 
strength (r) ^ which occur in vortex lattice representations, at the expense of 
more elaborate analysis. Three previously developed methods for improving the 
induced drag computation are those of Loth and Boyle (ref. 4), Goldhammer 
(ref. 5), and Clever (ref. 6). 

Loth and Boyle developed a theoretical Trefftz plane analysis of the 
wake of a single planform where the wake vortex sheet strength was assumed 
piecewise linear and continuous. However, the numerical results from 
reference 4 are seriously in error due to errors in the implementation of 
this method. More recently Goldhammer and Clever have developed near-field 
analysis codes that assume piecewise linear variation of singularity strengths 
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in the chordwise direction. Reference 5 assumed a piecewise constant span 
load, leading to a discrete vortex representation of the wake, while reference 
6 assumed a piecewise linear variation of span load, resulting in piecewise 
constant wake vortex sheet strengths. The current work utilizes the theo- 
retical model developed by Loth and Boyle, where the shed sheet strength is 
assumed piecewise linear, leading to quadratic variations of bound circula- 
tion and span load. 

A Trefftz plane analysis of a wing wake is performed using an assximed 
piecewise linear functional form of the wake strength, and solutions for the 
bound circulation for minimum induced drag are obtained using two optimi- 
zation techniques. First, Munk's criterion (ref. 7), applied at each wake 
segment midpoint, is used, similar to the technique developed in reference 4, 
but with important analytic differences in the handling of singularities at 
the ends of adjacent wake segments and at the wing tips, as well as the 
inclusion of variable wake segment spacing. Second, a direct optimization 
technique as discussed in reference 1 is developed for the assumed wake 
model. Results for both techniques are then compared with available exact 
solutions from Mangier (ref. 8) , Cone (ref. 9) , Lundry (ref. 10) , and 
Lundry and Lissaman (ref. 11) . 


SYMBOLS 


p,k 


matrix of coefficients in direct 
optimization technique 


A,B,D,E,F,G,J,K,R,T,U,W 

A' ,D' ,F' ,G' , J' ,K' 


constants appearing in normal wash 
expression (eq. 4) 


Ai . . , 


■ID 


’i: 


A 4 . . 


integrals appearing in normal wash 
expression (eq. 4) 


b 


wing span 


C 


L 


lift coefficient 


C induced drag coefficient 

D 


d 


vertical extent of circular arc 


D. . 
ID 


induced drag on segment i due to 
induced velocity of segment j 


G. 

1 


G. ,G. 


1 X 


variables containing unknown shed sheet 
strengths appearing in drag expression 
(eq. 7) 


2 



h. . 

ij 


h! . 

13 


’■^ 1 , j j '’■‘‘i, j 


distance between influenced point on 
segment i and influencing point on 
segment j 

distance between influenced point on 
segment i and influencing point on image 
of segment j 

integrals appearing in drag coefficient 
expression (eq. 7) 


k 

KirK2,K3 


I 


T 
R. . 

13 

r: . 
13 


s 


ref 


4 


s 

U 


w 


w 

n 


w* 



w 

o 


y 


induced drag efficiency factor 

constants appearing in integrals on p. 16 
and p. 28 

vertical extent of endplate or fence 

number of wake segments on one half of wing 

factors in normal induced velocity expressions 

projection of distance h. . onto the plane 
of influenced segment i 

projection of distance hi . onto the plane 
of influenced segment i 

reference wing area 

local wake segment coordinate 
wake segment half width 
free stream velocity 
total induced velocity 
normal wash velocity 

normal wash velocity induced at point 4=4. 
on segment i due to segment j 

normal wash velocity induced at point 4=4^ 
on segment i due to image of segment j 

constant appearing in Munk ' s criterion 
normal wash velocity expression (eq. 8) 

spanwise coordinate 


3 


z 


vertical coordinate 


^ b/2 

Y 

r 


r 

n 

A 

p 

o . . 

x: 


parameter for circular arc spanwise camber 

trailing vortex sheet strength 

bound circulation along the planform 
perimeter 

average bound circulation 
nondimensional spanwise coordinate 
Lagrange multiplier 
fluid density 


angle between direction of h . . and 

i: 

y-axis 


dihedral angle 


Subscripts 

i 

j 

n 

P 

o 


influenced point 
influencing point 
normal component 
diimmy index 
value at y = 0 


Superscripts 

' image segment quantity 

average quantity, or body axis 
coordinate 
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THEORETICAL DEVELOPMENT 


To effect a solution for minimum induced drag using the assumed wake model, 
the standard expression for induced drag evaluated in the Trefftz plane, given by 


C 


D 



I 

U U 


1 

cos <fl 


dy 


can be used. To develop expressions for the individual quantities in the 
integrand, i. e. bound circulation and induced normal velocity in terms of 
the wake vortex sheet strength, consider a view from downstream of the 
undistorted wing wake as shown in figure 1. The wake is broken up into 
N^ linear segments which can be of varying size. Uniform and cosine 
segment spacing are investigated in the current study. From the law of 
Biot-Savart, the total induced velocity at a point on wake segment i, located 
a distance -4^ from the center of that segment, due to an increment of 
shed vorticity from segment j, located a distance 4^ from the center 
of that segment, is equal to 


dw (-4 . , 4 . ) 
1 : 


Y (4 . ) d4 . 

2 1 

2iTh. . 
i: 


where 



+ (z , . - 4 . 
ID D 


and 



Q 

COS d) . + 4 . cos (b . ) 

D 1 1 

sin d) . + 4 . sin d) . ) / 

11 1 


z . , = z . - z . 
ID 1 D 
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Refer to figure 1 for details of geometry. The normal induced velocity is 


dw* (4 . . ) = dw (6 . ,6 . ) cos ((j> . - a, . ) 


n 1 3 13 

= dw (4 . f -6 . ) - — — 

1 : h. . 

1 : 




v/here a is the dihedral angle between points located distances 4 . and 
i 3 1 

4 . from the centers of segments i and j . is the projection of h^^ 

in the plane of segment i such that 


R. . = sin d) . (z . . - 4 . sin d) , + 4 . sin d) . ) 

i: IJ 3 3 1 1 


+ cos d) . (y . . - 4 . cos d> . + 4 . cos d) . ) 

1 13 D D 1 1 


Thus, the normal wash at 4^ on segment i due to segment j is 

, Y (A . ) \ cos A. (y. . - 6. cos A . + 4 . cos A. ) + sin (f>. (z. . - 4 . sin A . + 4 , sin A . ) > d4 . 

w* . (A . ) = ^ / 3 \ i ij j 3 i li! 1 i3 3 J i i f 3 

^ J (y, . - 4 . cos ^ . + 4 , cos ifi. ) ^ + (z. . - 4 . sin ^ . + 4 . sin ^ 

*'-S. i3 3 j 1 i 13 D D 1 1 


( 1 ) 


Next, consider the noirmal wash at a distance 4. from the center of segment 

1 

i due to the image of segment j, located on the right half of the wing. 
This normal wash, denoted by a prime, equals 


dw' (4 . ,4 . ) 
n 1 3 


-Y(4'.)d4'. cos(d). - a!.) 
3 3 1 J-3 


2tt 


hi . 
3-3 


-y(4!)d4'. R1 . 
■ 3 1 


2tt 


[h! .]2 
3-3 


where 


Rl . = sin d).(zl. -4. sind). +4. sind).) 
13 1 13 3 31 1 


+ cos d> . (y 1 . + 4 . cos d) . + 4 . cos d> . ) 
3- 3.3 3 33 . 1 


[hl.]^ = (yl. +4. cos d» . + 4 . cos d> ■ ) ^ 

13 13 3 3 1 1 

+ (z! . - 4. sin d> . +4. sin 6.)^ 
3-3 3 31 3. 
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The total induced normal wash at point i due to the image of segment j is 
then 


(2) 


w' (4.) 
n»3 1 



Y (4 . ) < cos (y! . + 4 . cos 4 . + 4 . cos 4 . ) + sin * . (z I . - i . 
J ( J- J 3 1 1 1 ID 3 

{yl . + 4 . cos + . + 4 . cos <^. ) ^ + (zl . - 4 . sin 4 , 

J-3 j j 1 1 13 3 3 


!in if. +4. sln4.)>d4. 
3 i f 3 


+ 4 . sin if . ) ■ 


where y ! . = y . + y . and z I . = z . . . The noimial wash at point i due to both 

ID 1 3 x: 13 

segment j and its image, equal to w* .(-6.) + w' (6 ), is then 

n ,3 X n,3 x 


w 


( 4 .) 


1 




ds . . 


The shed vorticity distribution, assumed piecewise linear as shown in 
sketch a, is defined for segment j as 


Y(4J = 


j+1 


+ Y. 


4 . 

-Jl 
s . 
3 




- Y 



V / 

Segment j 


Sketch (a) 


where is the value of the shed vorticity between segments j and j-1. 

The range of applicability of each equation for y is -s . to s . . Note 
that this assumption leads to piecewise quadratic distributions of bound 
circulation and span load. 

Then the normal wash at 4^ due to segment j and its image is 
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w . (-4 . ) 

n,3 1 


2 2tt 


r-"®jR. .di. 

/ ii"-J 

•'-s . 13 -s . 

3 D 


j R1 .dA . 

i3 1 


[ h : .]2 

■s. X3 


^ ^3+1 ' \ ±. 
2 2tt| 


p3 _ r ^"3 

*'-s. ij^j *^-s. 

I 3 3 


>6 .Rl .d& . 
.3 12 1 

[h' ] . 
i3 3 


( 3 ) 


Constants are defined as 


R., = A + B4. + 4. = A' + B4. 
13 3 1 3 


where A = y. . cos 4> . + z. . sin (J) 

13 1 13 1 

A' = A + 4^ 

B = -cos 4> . cos (f»^ - sin (})^ sin <{>^ = -cos (c|) ^ - <f>^) 

R'. = D + E4. + 4. = D' + E4. 

13 31 3 

where D = yl . cos 6, + zl . sin (f) . 

13 1 13 1 

D' = D + 4 . 


= cos (}) . cos (f)^ - sin sin (j)^ = cos(<j)j + <j)^) 


h?. = 4? + (F + 2B4.)4. + (G + 2A4 . + 4.) 
13 3 13 11 


= 4? + F'4 . + G’ 
3 3 


and 


where P 

F* 

G 

G' 


= -2 (y . . cos d) . + z . . sin d) ) 
13 3 13 j 

= F + 2B6 . 

1 

= z? . + y? . 

13 13 

= G + 2A4 . + 4? 

1 1 


(h'.)2 = 4^ + (J + 2E4.)4. + (K + 2D4. + 4^) 
13-1 11 11 - 


= 4? + J’4. + K’ 
3 3 
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where J 


= 2 (y I . cos 6 . - z I . sin d) . ) 

ID D 


J' = J + 2E4. 

1 

K = (z'j)2 + (y!.)2 


K' 


= K + 2D6 . + 4f 

1 X 


The normal wash is now written as 

^ ^ ^ —2 ^ *-lj> 


where 


Ai . . 

^iD 


A2. . 

^X3 




ID 


‘‘'^j A' + B4. 


di 


4? + F'A . + G' 
s. D 


-if 


’’’^j D' + E4_ 


-S . 

D 


4? + J’-A . + K‘ 
D D 


d4 


r'*‘®j A *4, + 

1 / 3 3_ 

/ 4? + F'4 . + G' 

-'-s. 3 3 


d6 


and 


A 4 . . = 
ID 


, D’4 . + £4? 

1 I _■ ■ D 3 

J + J'-* • K' 

m/ T T 


dA . 
D 


-s . 
D 


The above integrals contain constants, none of which contain 4^ 
be integrated analytically to give 


They 


1 j f .V . bf 


-)p + |(£n|s2 


? + F’s . + G' 
D D 


-* ln\ s^ - F’s . + G' I 

I 3 ^ I 


(4) 


can 
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where 


for 

If 4G 
Next, 

where 

for 

If 4K 
Also, 


P = 


^4G» _ f'2 


tan 


-1 


F' + 2s 


(4G* - 


- tan 


F' - 2s. 

-1| 1 

(4G' - F’2)1/2 


4G' - F'2 ^ 0. 


' - F' 


2 


0 , then P 



2 


F' + 2s . 

: 



4K’ - J'^ 7^ 0. 


- J-2 


0 , then Q 



2 

J' + 2s . * 
1 
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and 



Note that 

V4G' -f'^= |r+T4.| 

where 

R = 2 (y . . sin d) . - z . . cos i . ) 
ID D J-D D 

T - 2 sin ((|)^ “ (f>^) 

and 

^K* - J»2 = |u + W4^| 

where 


U = 2 (y I . sin d) . + z I . cos 6 , ) 

ID D ID D 

W = 2 sin (di , + d) . ) 

Now, written out in terms of the unprimed constants and the variable 
4^ , the integrals appearing in equation (4) for the normal wash are 


/ 



+ 2(A + Bs.)4. + 

D 1 1 


- - FSj + G 


+ 2 (A - Bs . )4 . + A? I 

3 X X 


-1 

(d - ^ + (1 - 

e2)^ ) 

^an ^ 

”j + 2s^ + 2E4^'j 

2‘rr 


U + W4 . 

X 


U + W4 . 

l- X 

J 


- tan 


-1 


J - 2s . + 2E4 . 
3 1 


U + W4. 

X 


- — £.njs^ + Js . + K 
4-rr ' j j 


+ 2 (D + Es.)4. + 4?| - - Js + K 

3 X X ' 3 J 


+ 2(D - Es.)4. + -4? 

D 1 X 


^ 3 - = 


^ II (BF^ - AF - 2BG) + (-F - 6AB + 4B^F)4^ + (4b3 - 4B)4? 


’ij 2 tts^ 


2 

I tan ^ 

> + 2s . + 2B 4. ' 
3. 1 

- tan“ 1 

"F - 2s . + 2B4 . " 
3 1 

M 




1 / 

|r + T4. 1 
1 1 

\ 

1 R + T4 . 1 
L ' 3 ' J 


|r + T4^| 

/) 


+ 




BF + (1 - 2B^)4^' 


llnls^ + FS. + G 
3 3 


+ 2 (A + BS . )4 . + 4^ t - ^*^1 s? - Fs . + 
3 11 3 3 


G + 2 (A - Bs . )4 . + 
3 1 


12 


csj -H 



A4. 


-1 

27rs 


(Ej2 - DJ - 2KE) + (-J - 6ED + + (4E^ - 4E)4? 




2 

|u + W4^| 



E 

TT 


1 


2tts . 
3 


D - EJ + (1 - 2e2)_4. 


Zn 


+ Js . 
3 


+ K 


+ 2 (D + Es . )4 . +6 
3 1 



Js. + K 


+ 2 (D - ES.)6 . 

3 1 



The four integrals above are now used to evaluate the induced drag on 
segment i due to the induced velocity at the wing due to segment j by 
integration of the product of bound circulation and normal wash. The induced 
drag for each linear segment is of the form 


D. . 

i: 



r(4.)(iw 

1 \2 n,3 X / 1 


and thus, the corresponding drag coefficient is 


C 


D,ij 


D. . 


pu-^s 


ref 



r(4. ) 

1 

u 


/ W . (4 . ) 
.P'3 L- 

\ “ 


d4 , 
1 


(5) 


The bound circulation is found by the integration of the shed sheet strength 
from the tip to the desired location, such as 


r(4. ) 

X 



Y (4 )d4 = r (-S . ) 
o X 



Y (4 . ) d4 . , 
X X 


That is , 


13 


r(4 


i> = rj-s.) . f (y.^^ + 3y.) . 


Y . , “ Y . \ -4 . 

1+1 ' X J _1 

2 / 2s. 

' 1 


where V (-s.) is the value of the bound circulation evaluated at 4. = -s., 

O 1 11 

Specifically 


i -1 r ^ / 

’ (-S.) = y] / = Y,s, + ^ Y|S + s \ + Y.s. ■ 

° ^ pti p’V ^ 


Then the induced drag coefficient on segment i due to segment j and its 
image is given by 


'« ^refA, I “ ""V \ 2U j i 


^i+1 - ^i\ 


2U / 2s^ 


Y + Y . 
1+1 

2U 


(A]^ . . + A 2 . . ) 
^ 1 : ^X3 




( 6 ) 


Define 


G. = 


r (-S . ) s . 
o 1,1 


+ TVT Y_. ,Y + 3y, 


i U 4U I ’ i+1 i 


Y + Y . Y . + Y . 

5 = 5. = J 


G. 

1 


2u : 


2U 


G. = 
1 


- ^1 : - Y 


2U 


G. = 
1 


2U 
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so that 


n (g.G.Ii. . + G.G.Io. . + G.G.Iq. . + G.G.I4. 

S^ef V 1 D 1 1 : 1 j 


+ G.G.Is. . + G.G.Ig . 

1 : ""x,] X J I 


where 


■ii,j = j ‘"'lij 


+ A 9 . . ) d4 . 

^3.J X 


l2 


psi 

i,j = / 

^ — Sj 


r+s^ 

^ 3 ^ -i “ / - (Al . . + A2 . .)<i6. 

i »3 J X ^XD ^13 i 

/ +Si 

-6 . ( Aq . . + Au . . ) d -4 . 

1 ID 1: X 

■^i 

(Aj . . + A2 , . ) 

42 ^d4. 

X 2 s. X 

•s- X 


l6 


(A 3 + a 4 ) 

= I -i? ^ ^d 4 . 

i/D / X 2s. i 

— s • 1 


( 7 ) 


These integrals with respect to 4^ are evaluated analytically using 
the MACSYMA symbolic manipulation language (ref. 12) as detailed in the 
Appendix. The following integrals require special consideration: 

1, Integrals of the type 

4^ £ki |42 + E4^ + D|d4^ n = 0,1, 2, 3 
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are evaluated exactly except where a logarithmic singularity occurs, i.e., 
at the endpoint. There, the singularity is omitted from the range of integra- 
tion. An analytical justification for the omission of such singularities has 
not yet been found. Instead, numerical studies have indicated that the omitted 
range can be varied from (10*"^) (s^) to (10“^) (s^) and still not affect the 
solutions for bound circulation of induced drag to four significant figures. 
This therefore provides some numerical justification for the restricted 
range of integration. 

2 . Integrals of the type 



-6^ 

1 


R + T^. 

1 


tan 


-1 


c + 2B6. 
1 

1 R + T 4 . 

' 1 


d 4 . 
1 


n = 0, 1 , 2 , 3 , 4 


are not generally integrable analytically; hence, they have been replaced by 
approximate integrals of the type 



+ K3) 

I R + T4 . 1 

' 1 


ci6. 

1 


which are then integrable analytically- The constants Ki , K2 , K3 , are 

chosen so the quadratic exactly equals the original inverse tangent function at 

4 . = ± s. and zero. When T = 0 and R 0 , the original integral is integrated 
1 1 

exactly. If both R and T equal zero, the remaining finite integral is 
of the form 



which can also be integrated analytically. For these three types of integrals, 
the denominator is checked for singular behavior. If singularities are found 
to occur, they are excluded from the range of integration. Again, numerical 
studies have shown that this does not affect the solutions for F or induced 
drag. Such singularities are found to occur only at the wake segment endpoints 
and only when considering the velocity induced by one wake segment on an 
adjacent segment. 
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OPTIMIZATION 


Using Munk's criterion (ref. 7) which states that 
w 

-r = w = constant, (8) 

cos 9 o ' ^ 

a system of linear equations is developed for the unknown shed sheet 

strengths, by equating the normal wash at each segment center to cos (p on 
that segment. That is, following the techniques of references 4 and 13, 



In the present formulation it has been assumed that the shed sheet strength on 

the wake centerline, equals zero. The resulting shed sheet strengths 

are scaled by the constant w . The ratio of w to the free stream speed is 

o o 

calculated from an evaluation of the lift coefficient. That is. 



Once the optimum shed sheet strengths are found, the drag is calculated 
directly as described in the previous section. A drag efficiency factor is 
found, defined as the ratio of (C ^/ttA) for a planar wing of equal span to 
the calculated induced drag for the nonplanar configuration. 

An alternative optimization technique which has been developed is to solve 
directly for the shed sheet strengths by writing the drag coefficient explicitly 
in terms of the unknown y's, as in the previous section, and finding expres- 
sions for the derivatives of with respect to the unknown scale factors. 

With 


T T 


: = 2 V V . . 

° i=i j=i 
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P'" 4 i=p-i ,^r!E ^(5 ,.+a.i2 1(5.13 ,. 

Jy) ®ref L 4 V 3 P-lO 3 P-1.3/ 2 \ j P-1,] 

®Mp 

-1.3) 


+ G.I 4 , . + G.I 5 - . + G.Ig 

J P-l/D 3 P-I/D : p- 


3 {C^ ... 




ref 


i^^P ^ ^P-1^ j ^ ^j^2p+m,j^ 


, m > 1 


The lift coefficient constraint is satisfied through use of a Lagrange 
multiplier in the objective function, so that the function to be extremized 
becomes 


N„, 

T T 


2 V y] c . + X ( 5^ c 
iti \p=i 


. — - c 

p U L 


Each equation in the optimization system can then be written as 


A , Y, + =0, p=l, . . . 

p,k'k L,p ' tr' > rj. 


where 


p,k S 


ref 


s + .s 
. P P -1 


+ I 2 


p+m 


,k-l^ 


N -P 
m=l \ 


^ p+m , k 


I 2 , + ll . 
p+m,k p+m,k-l 


2 ^^p-1 ^ ^P^ ^^^p,k ^ ^^p,k-l ^^p,k ^ ^^p,k“l^ 


4 r3p,k + ^‘+p,k-i ■ 
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' ^®Prk-l 




^^p-l,k 


^lp-l,k-l " ^^p-l/k 




l-pp-l.k ■" ^Vl.k-1 “ %-l,k %-l,k-l ^Vl.3^ 


- I - "v,p " ^V,P-1 - ^^lc,p - ^=k,p-l " - ^6 k,p-1 


4 i’^ 5 ^-l,P * ^^k-l.p-1 ■ ^“k-l,p 


^^k-l,p-l ^^k-l,p *'’ ^^k-l,p-l 


^^k-lrP ^^k-l,p-l 


k-1 


^^k-l,p ^ ^^k-l,p-l ^^k-l,p ^k-l,p-. 


"k^i ^i"k 


^^k,p ^^k,p-l ^^k,p ^^k,p-lj 


2 \ ^ ®k 


_1^ y ^k+m,p ^^k+m,p-l ^^k+m,p ^k+m,p-l^ 


And, since 


^L 


ref 


1 

J Z 

i=l 


) ■ St 
X X 


^i +1 


2 Y. 


i=l 


cos <i) . s . 

X X 


o 

U 
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then 



N^-P 

T 


ref 


E ( 

Tl=l \ 


COS (p s ] { s + s . \ + cos ^ s s 

^p+m p+m ) \ p p-1 ] ^p p p-1 


The (N^ + 1)^^ 


equation satisfies the lift coefficient constraint: 


N 

T 

E 

p=l 


C 


L,p' 


U 


0 


RESULTS AND DISCUSSIONS 


The theoretical development has been implemented in an FORTRAN computer 
program which is operational on a CDC Cyber 173 machine, but it is not 
currently documented or available to outside users. A measure of the conver- 
gence, accuracy, and effect of segment spacing for the two optimization 
procedures is shown in figure 2. The induced drag coefficient calculated 
from direct integration of the drag expression obtained from the shed sheet 
strength values is compared with the exact result for a planar wing. Results 
are shown in terms of an efficiency factor, k, which is defined by Cone 
(ref. 9) and Lundry (ref. 10) as 

C exact planar wing 
k = 

C calculated 
D 

It is seen (note the stretched scale) that for both optimization techniques, 
k rapidly approaches the desired value, 1.0, as the number of wake segments 
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increases. Induced drag efficiencies for the direct optimization technique 
are very close to those calculated using Munk's criterion. However, as seen 
in figure 3 for = 10, there are slight differences in the shed sheet 
strengths near the wing tip which are especially noticeable for equally spaced 
wake segments. Results similar to those shown in figure 3 indicate that 
the behavior of solutions for = 25 and 50 is similar to those for 

= 10, except the y values lie closer to the exact distribution. Note 
the dramatic improvement in the representation of the wake vorticity distri- 
bution that is achieved using cosine spacing of the wake segments. Since 
the two optimization techniques yield nearly identical results, the remaining 
comparisons of the current theory with previous work will display only the 
Munk criterion optimization results. 

In figure 4 induced drag efficiency factors for a planar wing computed 

using the current theory are compared with similar results obtained using a 

discrete trailing vortex Trefftz plane drag optimization program developed by 

J. R. Tulinius, B. B. Gloss, and J. L. Thomas of the NASA Langley Research 

Center using the method of reference 14. The error in the present theory 

with N^ = 10 is 1 percent with equal spacing and 0.2 percent with cosine 

spacing. This error using equally spaced segments is seen to be approximately 

one-fourth to one-fifth the error in efficiency obtained using the technique 

of reference 14. Errors for the current theory using cosine spacing are five 

times smaller still. The bound circulation values calculated from the Munk 

criterion solution shed sheet strengths are compared with the exact elliptical 

distribution in figure 5 for N =10 and 25, Errors in calculated T values 

T 

are even smaller than any errors in y. 

Optimization results for nonplanar configurations are compared with exact 
solutions in figures 6, 7, 8, and 9. The k values for a nonplanar configura- 
tion composed of a flat wing with a vertical endplate are compared in figure 
6 with those of Lundry and Lissaman (ref. 11) . Results shown are from the 
downwash criterion optimization using N^ = 30, with the wake segments 
distributed nearly uniformly on the wing and endplate and tj) = 89.9°, and with 
cosine spaced segments and (p = 89.7°. The endplate is only approximately 
vertical to avoid numerical difficulties at = 90°. The k values are 
somewhat sensitive to the value of (p chosen, but at <|> = 89.7° with 30 cosine 
spaced wake segments the accuracy of the present method is comparable to the 
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method of reference 14 using’ 100 equally spaced discrete vortices. Note the 
relatively large error at Z/s - 0.2 for equal spacing. 

In figure 7 the induced drag efficiency for a wing with constant nonzero 
dihedral (<j) = 30®) outboard of ri = 0.5 is shown compared with the solution 
given by Lundry (ref. 10) at Z/s = 0. Results also compare well with a 
solution given by Mangier (ref. 8). The current code cannot handle vertical 
fences (Z/s >0), but this capability could be added. 

In figures 8 and 9 results for wings having spanwise camber such that 
the wake is an arc of a circle are compared with an exact solution for such 
wakes by Cone (ref. 9) . All results are from the constant downwash formula- 
tion of the optimization. Figure 8 compares the induced drag efficiency, and 
figure 9 shows the bound circulation values. The parameter B is defined as 


6 = 


d 

b/2 


where d is the maximum vertical dimension due to the spanwise camber. A 
value of 3=0 corresponds to a flat wing, and 3 = 1-0 corresponds to a semi- 
circular wing. Present results with N^ = 25 again display an accuracy 
comparable to that of the method of reference 14 with 100 discrete wake 
trailing vortices. 


CONCLUSIONS 


The current drag optimization techniques assuming a piecewise linearly 
varying wake vortex sheet strength agree well with each other and with 
available exact solutions for a variety of Trefftz plane wake shapes. 
Agreement for overall drag efficiency factors and bound circulation distri- 
butions is generally better than one percent for on the order of 25 to 50 
wake segments. Computational times generally run from 1 to 10 seconds on 
a CDC Cyber 173 machine. Accuracy of the present theory is approximately 
four to five times better than a discrete trailing vortex theory using the 
same number of unknowns and equal spacing of the wake segments. Cosine 
wake segment spacing leads to a further increase in accuracy. 
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APPENDIX 


EVALUATION OF INTEGRALS IN DRAG EXPRESSION 

The following is a compilation of the integrals utilized in the 
evaluation of the drag expression [eq. (6)] in the text, which are not generally 
available in integral tables. They have been evaluated through use of the 
MACSYMA symbolic manipulation language (ref. 12 ) , and are repeated here for 
completeness . 


^ + E-6 + D)d6 depend upon 

-s 

whether E^ - 4D is positive, negative or zero. Results for all three are 
given, starting with E^ - 4D > 0. For this case 


^ + E4 + D)d6 = ”jVe^ - 4d In^ 


2s - V E^ - 4D + E 
2s + V E^ - 4D + E 


+ s Zn (s^ + Es + D) (s^ - Es + D) - 4s 


E s^ + Es + D 1 “ p I -2s - \E^ - 4D 

+ “ • tn + — V - 4D I - 7- — 

s^ - Es + D \-2s +Ve2 - 4D 
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For 


< 

f 

-s 


InU^ + E4 + D)d6 = — 

o 


1 E^ - SDE-^ + 4D" 


V 


E^ - 4D 


- tni 


'2s 


Ve^ - 


4D + E 


■V 


2s +VE^ - 4D + E; 


+ Znl 


-V 


-2s - \E^ - 4D + E 


\-2s 


+ yjs^ - 4D + E/ ) 
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- 4D 
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4D + E 


4D + E, 
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+ In 


( S^ - Es + d \ 
s^ + Es + D/ 


+ ^ fa p + ES d \ ^ E|1 
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(3D + E^)s E 
2 


- 4D = 0, 
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InU^ + E4 + D)d6 = S In (s^ + Es + D) (s^ - Es + D) 


+ E In 


( E + 2s \ 
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4s 
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n . 0 V /> /s^ + Es + D \ E^ - 2 d 
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tnu^ + E4 + D)d4 = s-tn(s^ + Es + D) (s^ - Es + D) 


Zn f- + P \ _ 4s _ -W4D - • 
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r" ^3 ^ D)d4 = ? gi 

J -S 



- Es + D 
+ Es + D 



-h Es + D 
- Es + D 


+ 


- E^ + 6DE^ - 8D^E 


V4D - 




_ tan-l/ ^ - ^(3D - E^) 

6 2 

/ 4 D - E^' " 


2 . Integrals of the type 


s 



integrals 




are in general replaced by approximate 


+ K2-6 + K3) 
R + T 4 


di , where the 


Kl, 


K2 , K3 are chosen to force 


the quadratic approximation to the inverse tangent function through the exact 
integrand at s = -s, 0 , +s. Hence , the integrals are evaluated as follows: 

For n = o. 



+ K2^ + K3 

I R + T 4 1 


d6 


(K3T^ - K2RT + KiR^) 


Zn 


;R + TS 


R - TS 


+ — {2K2T - 2 KiR) 
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for n = 1 , 


o' -s 


S + K2-6 + Kg) (KsRT'^ - K 2 R^T + K^R^) /|r - Ts 

4 d 6 = Zn I ^ — 


R + T4 


R + TS 


+ — — (6K3T^ - 6K2RT + 6KjR^) + 
3T^ 


2Kj 

3T 


for n = 2 , 
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In 
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6 t 2 t 4 
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for n = 3 , 
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for n = 4 


/ 


s (K 6 ^ + K 26 + K 3 ) (KsR'+T^ - K 2 R^T + KiR^) 


d 4 = 


-s R + TA 


R + TsI \ (KoT - KiR) 

, ln{ 1+2 ^5 
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When T = 0 and R ^ 0, the integrals are evaluated analytically at 


n = 0 , 
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continued 



When both R and T = 
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Figure 2. Effect of optimizing procedure and spacing 
on convergence of the induced drag 
efficiency for planar wings. 





spacing on the wake strength 





Figure 4. Convergence of induced drag optimization 
using Munk's criterion for planar wings. 











